Mean-Field Nematic— Smectic-A Transition in a Random Polymer Network 



Peter D. Olmsted* and Eugene M. Terentjev 
Cavendish Laboratory, Madingley Road, Cambridge, CB3 OHE, U.K. 
(February 1, 2008) 

Liquid crystal elastomers present a rich combination of effects associated with orientational 
symmetry breaking and the underlying rubber elasticity. In this work we focus on the effect of 
the network on the nematic-smectic-A transition, exploring the additional translational symmetry 
breaking in these elastomers. We incorporate the crosslinks as a random field in a microscopic 
picture, thus expressing the degree to which the smectic order is locally frozen with respect to the 
network. We predict a modification of the NA transition, notably that it can be treated at the 
mean-field level (type-I system), due to the coupling with elastic degrees of freedom. There is a 
shift in the transition temperature Tna, a suppression of the Halperin-Lubensky-Ma (HLM) effect 
(thus recovering the mean-field continuous transition to the smectic state), and a new tri-critical 
point, depending on the conditions of network formation. When the nematic phase possesses 'soft 
elasticity', the NA transition becomes of first order due to the coupling with soft phonons in the 
network. We also discuss the microscopic origin of phenomenological long-wavelength coupling 
between smectic phase and elastic strain. 



I. INTRODUCTION 

Randomly crosslinked networks of polymer liquid crystal (LCP) materials (liquid crystalline elastomers and gels) 
have been a subject of substantial experimental and theoretical activity in recent years. A newcomer to this area can 
find out about about the synthesis of side-chais. and main-chain systems, characteristic physical effects, and concepts 
related to the nematic state in review articleta-EI. Three major factors determine the behavior of these remarkable 
materials; liquid crystalline symmetry breaking; rubber elasticity coupled to the resulting anisotropy, producing a 
highly mobile principal axis; and, finally, randomly placed (and sometimes randomly oriented) network crosslinks. 
Random disorder introduced by these crosslinks encourages elastomers formed in the isotropic state to cool into 
polydomains, i.e. highly non-uniform textures of their liquid crystalline phases. If the network is formed in the 
presence of an external field, or in a uniform monodomain low-temperature phase, the resulting elastomer remains 
macroscopically uniform. ^ ^ 

Smectic elastomers (refs.Elu represent a few recent examples) comprise liquid crystalline polymers crosslinked into 
a network which contains liquid crystalline groups ordered into a one-dimensional density wave, or a smectic state. 
One may consider either main chain or, more commonly, side-chain LCP's. The smectic order is typically obtained 
by forming the network in a smectic state {e.g. by crosslinking at temperatures below the nematic-smectic-A, or NA, 
transition temperature in the melt). It can also be obtained by cooling a sample that has been prepared in the 
nematic state through the transition temperature Tjva, where we expect Tjva < T^a- The NA transition in networks is 
expected to be qualitatively very different from the same transition in melts because of the presence of elastic strain 
degrees of freedom and the local effects of random crosslinks. 

In the conventional NA transition, the degrees of freedom are the smectic order parameter V'(r) and the director 
fluctuations dn. The parameter ?/'(r) describes the departure-jof the mesogen center-of-mass density p(r) from a 
uniform density po, in a form of single wavelength modulationQu 

p(r) = po |l + -L [^(Oe^-i-- + c.c] I . (1.1) 

The nematic state is parameterized by the director n, which indicates the axis of preferential alignment of the 
mesogenic groups, and 5n indicates the fluctuations of this axis. While the NA transition should be continuous 
according to symmetry arguments, director fluctuations may change this picture. Halperin, Lubensky, and MatI 
argued that the coupling to the fluctuating "gauge" field 5n should induce a first order transition in type-I smectics 
(and the analogous superconducting system). In type-I smectics the characteristic length for penetration of the director 
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twist into the ordered smectic state is much smaher that the correlation length for the smectic order parameter, so 
that director fluctuations may be treated at a mean field level. Type-II smectics and superconductors, in which the 
gauge field fluctuations are correlated on length scales comparable or larger than that of the order parameter ip{r), 
require a more sophisticated anal^^is. The beliefs about the nature of this transition have varied over the years. A 
renormalization group e-expansionfl yielded no fixed point for physical order parameters, suggesting that th£j transition 
is weakly first order. However, duality arguments combined with Monte Carlo calculations on a latticecl predicted 
a continuous transition in the universality class of the inverted 3D XY-transition (that is, the amplitude ratios -are 
inverted) , experiments have generally yielded a continuous transition of the non- inverted 3D XY classEj E3 (but seellj) . 
Tonenlil has_argued, on the basis of a dislocation-melting theory, that the transition should be continuous; Andereck 
and Pattonli3 have calculated critical exponents within self-consistent perturbation theory; and Radzihovskyllj has 
performed a similar calculation for the aaalogous type-II superconducting transition, demonstrating the existence of 
a fixed point not found in e-expansion ofl. 

While our understanding of the NA transition could be incomplete, the nature of the smectic state is well- understood. 
The smectic one-dimensional density wave, of wave number qo, suffers the Landau- Peierls instabiliiji-and the system 
exhibits quasi long ranged order (with an algebraic decay of density correlations at long distances )EIl'ta. The source of 
this instability is the Goldstone modes corresponding to the arbitrary phase of the smectic density wave, whose long 
wavelength free energy is given by the Landau-Peierls elastic energy of a one-dimensional solidEj. In three dimensions 
this energy yields algebraic decay of correlations in the ordered state, analogous to systems with broken continuous 
symmetry in two dimensions. 

In contrast to a liquid crystalline melt, a liquid crystalline elastomer has an additional set of long-wavelength degrees 
of freedom: the elastic deformation v(r), defined by the local network displacement, r — > r -I- v(r). Smectic order 
couples to the elasticity through the crosslinks, and this coupling is manifested in two effects: 

A. Crosslinks pin the smectic phase and break translational invariance. The mesogens are either incorporated into 
the polymer backbone (main chain LCP's) or tethered to the network (side chain LCP's), and the stretching 
energy for relative translations between elastic displacements of the crosslinks Vz (r) along thi5 direction of the 
smectic density wave and the layer displacement u(r) can be written phenomenologically asEB 

2Fd = A [d\[vz{r)-u{v)]\ (1.2) 



This energy is defined on lengthscales of order or longer than the characteristic mesh size and, most importantly, 
is present for uniform relative translations. This coupling restores the Bragg peaks associated with smectic 
ordered. In Section II we argue that this terrpJs in fact only a metastable term, and the relative displacement 
^^(r) — u{r) may be relaxed by layer hoppin^J with a characteristic relaxation time. 

B. The preferential reduction of mesogen density around a crosslink due to, e.g., the steric exclusion of the mesogen, 
leads us (see Section II) to model the crosslinks by a local random field which adjusts the smectic phase: 

F„p = 7 [d\cir)\i^ir)\cos[qo{z-u{r) + vz{r))], (1.3) 



where c(r) is the local crosslink concentration. 

The smectic elastomer thus contains four contributions to the continuum free energy in addition to the smectic and 
nematic free energies associated with the NA transition in liquids. These are: 

(i) The elastic free energy of a uniaxial solidH, written in terms of the symmetric elastic strain eap = ^(Vaf/j-t-V/jWc). 
This contribution describes the phonon field (in most elastomers - incompressible) in the rubbery network. It 
couples to the relevant variables of the NA transition through the following effects. 



(ii) Eq. (1.2), the term penalizing relative shifts in the smectic phase variable u(r) and the phonon displacement 

parallel to it, Wz(r). 

(iii) Rubber- nematic couplingsil, penalizing relative rotations of the elastic strain 11 = i(V x v) and the nematic 

director, u) = n x Sn. 



(iv) The random field term, Eq. (1.3), describing the effect of crosslinks on the phase of the smectic order parameter. 
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In this jwprk we consider the microscopic origin of Eq. (|l.2|) and treat the random field in Eq ( |l.3| ) by the repUca 
formahsnO, with the following results [embodied in Eqs. (^:.6| )] for the mean-field theory of the NA transition in the 
network: 



1. Near the putative continuous NA transition Eq. (L2) may be ignored, since it is primarily a dynamic effect. At 

low enough temperatures (when the characteristic timescale becomes essentially infinite) this term must be 
considered. 

2. At the mean field level ncmatic director fluctuations no longer induce the Halperin-Lubensky-Ma first order smectic 

phase transition: they acquire a mass which reduces their effect on the NA transition. The existence of this 
mass in fact makes the mean field treatment (and the associated type-I assumption) essentially exact. 



3. However, in the special case where the nematic elastomer possesses soft elasticitM-Lsee Eq. B.5 below), characteristic 

of spontaneously breaking the symmetry from the isotropic to nematic statesE3E3, elastic strain fluctuations [the 
phonons v(r)] restore the HLM effect and the concomitant first order phase transition. 

4. If the nematic state is field-induced rather than spontaneous {i.e. the network has been formed to record the 

broken orientational symmetry) , the corresponding network phonon modes are conventional and the type-I NA 
transition should be continuous. We predict a tri-critical point in the crossover region between these two regimes. 

5. The effect of disorder on the NA transition is a simple renormalization of the transition temperature Tjva- The 

transition temperature increases because crosslinks localize the smectic phase variable and encourage order in 
the disordered phase. However, slightly below Tjva the smectic state crosses over to a 'glassy' state characterized 
by replica-symmetry-breaking. 

6. We emphasize that the typical nematic or smectic network is most likely to be 'hard', and that the 'soft' case 

we discuss extensively here is, for most practical cases, a theoretical construct which provides a framework for 
understanding the implications of a smectic state coupled to an underlying elastic continuum. 

Our description of the NA transition in the network is mesoscopic because we treat the smectic and director degrees 
of freedom as smoothly varying fields while including the microscopic interaction of the crosslinks with the smectic 
order parameter. The calculation in this paper sets up the proper coarse-grained model with which to investigate the 
effects of disorder on the smectic network, which we discuss briefly here and in more detail in a forthcoming work. The 
coarse-grained nematic and smectic fields are defined on a length £ic of order a correlation length for liquid-crystalline 
fluctuations, while the elastic strain field is coarse-grained on lengths down to the mesh size £x in the network, which 
may be larger than £ic. These different lengthscales must be kept in mind throughout. 

This paper is organized as follows: in Section II we present a more detailed discussion of the two primary effects of 
the network-smectic couplings; in Section HI we present the model free energy for the system; in Section IV we sketch 
the integration of the elastic strain and director fluctuations and present an effective free energy for the smectic in 
the presence of disorder; and we conclude in Section V. Our primary results are the effective Landau energies at the 
end of Section IV and the discussion in Section II about the smectic-elastic couplings. Our goal here is to understand 
the various terms in the Hamiltonian describing the NA transition, and identify the effects of disorder and rubber 
elasticity on the instability to the smectic state. We leave for the future a detailed analysis of the effects of disorder 
on the low temperature smectic state. The appendix contain some technical details of the replica calculations. 



II. SMECTIC-ELASTIC COUPLINGS 



A. Random Field 



In this section we justify the random field coupling, Eq. (1.3). This interaction represents the pinning of the smectic 



phase to network inhomogeneities, which we represent by a random distribution of crosslinks. Let us focus on side- 
chain LCP's, Fig. 1, since these constitute the most commonly synthesized liquid crystalline elastomers. The three 
constituents of the elastomer are flexible backbone monomers, crosslink groups, and mesogenic side-groups attached 
to the backbone with flexible spacers. In elastomers (as opposed to tough and brittle densely-crosslinked resins) the 
volume fraction of crosslinks is much smaller than that of the backbone and mesogenic groups, typically less than a 
percent of the aggregate. Hence we treat the environment of the mesogens as a uniform mesogen/backbone mixture 
with interspersed crosslink groups. 
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FIG. 1. Sketch of mesogenic groups in a smectic elastomer. 



It is reasonable to expect a steric repulsion between the mesogen and the crosslink which is enhanced relative to the 
steric repulsion between the backbone and the mesogen, simply because four chains come together at the crosslink. 
The presence of a Flory-x parameter between the mesogen and backbone groups contributes a similar effect. While 
there will also be an interaction with the local nematic order, since the mesogens may adopt an alignment dictated 
by the crosslink, we concentrate in this work on the effects of disorder on a smectic phase appearing out of a uniform 
conventional nematic state. Hence, we model the crosslinks by a local random field which induces smectic order and 



fixes the smectic phase, leading to Eq. (1.3). The corresponding energy contribution for each crosslink in the smectic 
potential is 

Ffl^ =^7l^(R«)|cos[(Zo(2»-u(R»))] , (2.1) 

i 

where is the position of the j*^ crosslink. The coupling constant 7 can be estimated on the basis of Fig. 1. We 
shall assume that the crosslinking point localizes the positions of <j) monomers (with (p the crosslink functionality). 
The barrier for such an object to 'tunnel' through a smectic smectic layer is a molecular characteristic of the material 
describing the degree of miscibility of backbone and the mesogenic side-groups, and is roughly of order (px, where x 
is the the Flory-x parameter between mesogens and backbone, which includes both steric and energetic interactions. 
Introducing the continuum crosslink density c(r) = y^ „- ^(r — Ri), so that under an elastic distortion R^ — s- R^ + v 



they distort by c(r) — > c(r — v), we can transform (2^) to collective variables. After changing variables r' — r — v(r) 
we obtain 

F^F^l /Ac(r)|V'(r + v(r))|cos[qo(2-u(r) + w^(r))]. (2.2) 



(Changing variables in the argument of u(r) introduces higher-order gradient corrections which are irrelevant in a 
mean- field treatment). 

Rather than working with the_discrete crosslink positions, we represent the crosslink concentration c(r) by a 
Gaussian probability distributionEZi: 

P[c] ocexp(-i /dV^l , (2.3) 



,27 2N^ ^ 

where Nx is number of crosslinks per unit volume of the system. For this random distribution of crosslinks the 
characteristic moments are (c(r)) = and (c(r)c(r')) = N^Sijc — r'). 

One should note, at this point, that a different situation emerges for a network formed in a deep uniform smectic 
phase: the polymer backbone in such a system is highly constrained between the smectic layers and, therefore, the 
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crosslink distribution (2.3) has a positionally modulated kernel. This possibility will have a profound effect on the 
phase ordering in the smectic phase, but is not relevant for the NA trpHsition description we are concerned here. Also, 
in a physical system crosslinks fluctuate about their mean positiong2a, which broadens the kernel of Eq.(2.3) into a 
Gaussian distribution with a width (in real space) proportional to the extent of the typical crosslink fluctuation. We 
ignore such fluctuation effects for now. 



B. Translational Invariance 



The second effect of the network on the smectic phase is due to stretching the polymer backbone. This lea ds t o-a 
free energy cost for uniformly displacing the smectic phase relative to the elastic displacement, given by Eq.(1.2)Ej. 
Whereas the contribution of the previous sub-section acts only at the crosslink position, the stretching effect influences 
the smectic phase all along the strand between crosslinks, and is coarse-grained at the mesh size or larger. To 
understand this in detail we first note the distinction between the smectic phase variable qqu and the z-displacement 
rji of the center-of-mass of the i*'' smectic mesogen. 

The phase qou is a coarse-grained variable and the displacements {iji} are microscopic variables. To calculate the 
position of the smectic phase from a microscopic picture one must, in principle, average over all the positions of the 
mesogens. ia an equilibrium smectic liquid the mesogens lie in a smectic potential given by, for example, a cosine 
modulationtj, but are not fixed in one trough of this potential: rather, they fluctuate back and forth over the barriers. 
In a strongly-ordered smectic state the activation over the barriers is very rare, while in a weak smectic state these 
fluctuations are quite common. The 'motion' of the smectic phase thus corresponds to the average motion of the 
mesogens' centers-of-mass. 

Now consider an elastomer in a smectic state and imagine displacing the smectic phase while fixing the crosslink 
positions. For clarity we focus on a side-chain liquid crystalline network, but the argument applies to a main- 
chain network as well. The displacement of the smectic phase displaces the average positions {rji} of the mesogenic 
sidegroups. Since these are tethered to the polymer backbone, this costs roughly the entropy of displacing the 
center-of-mass of a polymer chain a distance u while keeping the endpoints {i.e. crosslink positions) fixed. A simple 
calculation leads to an energy per chain of Vu ^ ksTv? / {No?), with N the number of monomers of size,-* between 
crosslinks. In addition to this elastic force, an individual mesogen feels a force due to the smectic potentialE3. At zero 
temperature this potential enforces the separation of the smectic phase variable and the crosslink displacement, but 
at finite temperature this barrier may be overcome by hopping. By hopping we mean that the mesogen fiips from one 
smectic trough to another, in such a way as to relax the strain. This allows the smectic phase variable to increase 
without bound while roughly localizing the mesogens' centers-of-mass. 

Hence a timescale r* separates liquid-like smectic behavior at long times from 'solid-like' smectic behavior at short 
times, the latter characterized by the elastic modulus A. For temperatures above the nematic-smectic transition 
temperature r* is extremely small and we may ignore this elastic effect. However, as the network is cooled deep into 
the smectic state r* grows, and for all practical purposes we must include this term. Since we are only concerned 
here with very weak smectic phases in this work, we leave further discussion of the interesting dynamics of a smectic 
network to another work. 



III. THE MODEL 



In this section we introduce the full free energy of the system, and set up the replica calculation. The NA transition 
in a liquid is described by the Landau-de Gennes free energy, 

2F„,, = jd\ {T|^(r)|2-Ki/3|V'(r)|4+gj_|(Vx-zgo<5n)^(r)|2-K5|||9,^(r)p}, (3.1) 

where the minimal coupling (V_l — iqaSfi) satisfies rotational invariance of the director and the layered systcmB. In 
the high-temperature nematic phase director distortions Sh are penalized by the Frank free energy, 

2Ffr = Jd^r {Ki{V-Shf + K2{h-VxShf + KsihxVxShf} . (3.2) 

In this work we use the one constant appr|ejtimation, Ki = K2 = = K . The elastic energy of an underlying 
uniaxial solid can be written in many waysE3, and we choose the representation through the traceless strain tensor 
iap = £0/3 — 5Tr[e](5Q./3, in order to deal explicitly with the case of a nearly incompressible network: 
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j d\ [al, + 2C2e-..Tr[e] + C3 (Tr[6])' + 2^4(6-2, + 2el^ + e^^) + AC^iel, + j} , (3.3) 



where the z-axis is chosen along the nematic director. In the isotropic case, the elastic moduli {Ci} transform to the 
Lame coefficients as Ci — 2fi ; C2 — 0; C3 = A + |/i and C4 = C5 = /i. In most cases we take the limit C3 00, 
since A//i ~ 10'' — 10^ in a typical rubber. The rubber shear moduli scale as /i ^ N^hgT, where ~ is the 
number of crosslinks per unit volume. 

The rubber-nematic freei-ejiergy penalizing relative rotations of the director u) ~ h x Sn and elastic strain ft = 



^(V X v) can be written asL— 

2F„„„ = Jd^- ^^bl[{n-(v)xhf + 2b2n-e-{n-uj)xhj , (3.4) 

where 61 and 62 are proportional to, respectively, the nematic order parameter and Q, multiplied by N^ksT. 
A positive coupling 62 corresponds to a network which favors parallel alignment between mesogen and backbone 
orientations, while a negative &2 favors perpendicular alignment. 

For the special case of nematic elastomers which possess an arbitrary choice of the nematic axis (because of a 
spontaneously broken symmetry of a network, formed in the isotropic state) rather than a 'quenched' axis of alignment 
dictated by, for exampley-|aa applied field during crosslinking or crosslinking in the nematic state, the relationship 
C5 — &|/(46i) = holdaSaO and the phonons described by ixz and iyz are "soft" {i.e. the corresponding phonons 
have fluctuation spectrum {\u{q)\'^) ^ g^* instead of the conventional behavior). The deviation of the network 
from this curious soft case is thus parameterized by 

(3.5) 

The presence or absence of these modes plays an essential role in describing the NA transition. While a perfect soft 
nematic elastomer (A = 0) is unlikely to be found in typical experiments due to random stresses in the nematic state 
or a predetermined alignment direction during network formation, wc include the soft case to present a complete 
theoretical picture. 

The partition function for the system is then given by 

Z = y V^/jVvVSh exp-^{F^^ + + F,i + F„_, + F^^]. (3.6) 

To find the effects of crosslinks on the NA transition in the network and, eventually, on the nature of the smectic 
state, we average over the disorder associated with these crosslinks. To compui|e-|quantities such as the quenched 
disorder-averaged free energy and correlation functions we use the replica trickcj to write the free energy of the 
system as 

F/k^T = -(logZ)p = lim ^ ^, (3.7) 



introducing n replicas of the system. The disorder average over the distribution P[c\, Eq. (2.3), couples the replicas 
together into the term FrepV- 

2 AT " f 

2Fr,,i ^ J2 d'r\r{r + ^r^{r))\\^\r + v\r))\ (3.8) 

^ a,b=l'' 

X cos [go {z - M°(r) + < (r))] cos [qo {z - u^{r) + vl{r))] . 

Note that u fluctuations are implicitly included in the measure for ip fluctuations. From here on we consider a 
uniform ■0 for the purposes of understanding the mean-field behavior of the transition, and as such will not consider 
fiuctuations of its phase. 

We next coarse-grain the system by averaging over the period of the smectic modulation qQ^ along the layer normal 

— 1/2 

z. This procedure requires that both the 'mean-free-path' between crosslinks £x ^ Nx and typical length scale of 
the relative translations u{r) — Vzi'r) be large compared to the layer spacing 2TTq^^. We find 

n „ 

Frepi = -r ^ / d\ irW^'l COS [<zo («'^(r) - v\r))] , (3.9) 

a, 6=1 •' 
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where F = T^iVa / (2kBT ) . T he free energy of the system is now given by Eq. ( |3.7| ), with n rephcated copies of the 
partition function Eq. ( |3.6[ ), with F^jp replaced by F^epi- Note that the similar form of the replica Hamiltonian, 
containing a translationallu-jinvariant cosine of fluctuating field {v°' — v^), appears in the problem of random flux 
pinning in superconductorSLj. 



IV. MEAN FIELD PHASE BEHAVIOR 



To see the effect of the crosslinks on the NA transition in the network we integrate out the strain and director 
fluctuations and examine the stability of the resulting effective theory for a uniform \ip\. This procedure is valid 
when the penetration length A for director twist into the smectic is much smaller than the smectic coherence Jength 
^ = (ffJ-/''')^^^ (tJ^is situation is often referred to as the type-I smectic-A, by analogy with superconductors^. For 
an elastic network, nematic fluctuations have a "mass" m « <ZoffJ-|V'P + biA/C^, where the first term is the mass of 
due to smectic fluctuations and the second term, ~ N^hgT, is due to the coupling to the clastic network. As the 
transition is approached we have ^ P/t and the twist penetration length is given by A = {K/m)^/^, so that 



A 



tK 



1/2 



(4.1) 



This Ginzburg parameter vanishes as the critical point is approached (t — s- 0), so that the condition for a type-I 
smectic (A/^ <C 1) is trivially satisfied and the mean- field treatment here is exact. [Note that for the perfect soft 
system, A = 0, the system may or may not be type-II, as with ordinary smectics.] Elastic fluctuations enter the 
theory at Gaussian level and may be integrated out at a mean-field level. A self-consistent check on the neglect of 
these fluctuations is to ensure the smectic correlation length at the fluctuation-induced first order transition is of 
order or larger than both the mesh size (to satisfy the coarse-grainijBg procedure) and the twist penetration length. 
To perform the calculations we choose a convenient geometry (cf.Q). We decompose the displacement into 



(4.2) 



where Vz is parallel to the nematic director; v± is normal to the director and belongs to the plane defined by the 
director and the wavevector q in Fourier space; and Vt is in the direction defined by q x n. Similarly, Sh = {6ht,5h±). 
Integrating out the director fluctuations is straightforward since they are not involved in the random coupling, and 
we find 



(4.3) 



{Z^)P = ff[ VrVw'^ exp --^ I ^ [Fi + Fel (V^ + Prepl 



The effective smectic free energy density is renormalized to 



V 



irlV-P + im'' + ksT / log h + Doiij) + Kq 



(4.4) 



where V is the system volume and we have used the one-constant approximation for the bare Frank elastic constants 
Here and below, / = J (Pq/ {2tt)^ . The parameter Dq — 
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_qQ is proportional to the square of the smectic order 



parameter, and the logarithm in Eq. (4^) is the determinant of the quadratic form from integrating the director 
fluctuations. The factor Dq{4>) in this term gives rise to the Halperin-Lubensky-Ma effect in a mean- field type-I 
system. That is, when the director is not coupled to the network (6i — 0), director fluctuations are massless and 
the g-integration rcnormalizes the transition temperature T^^, and yields a negative cubic term which induces a first 
order transitionB. However, massive director fluctuations (6i ^ 0) imply only even powers of l^/^l in the expansion of 
the integral, destroying the HLM effect. 

Eliminating Sn also rcnormalizes the rubber elasticity, yielding Fei, which is a quadratic form in the phonon 
variable v. Together with F^epi, the effective energy governing the elastic degrees of freedom is given by 



n 



J v''(q)-M(g,V').v'^(-q)-F^ J d^rW^'lcosqoiv: - v",), 



(4.5) 



where the kernel matrix M^q, tp) in F^i is given in the Appendix. 
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The next step is to integrate out the phonon variable v(r) and extract the effective free energy as a function only 
of tjj. Since Ti. is non-Gaussian this is a non-trivial step. Wer^rform this integration by applying the Gaussian 
Variational Method (GVM) as introduced by Mezard and ParisiKl for random systems with translationally invariant 
(in replica space) Hamiltonians, like our cos[w° — v^]. In performing this integration we assume replica symmetry for 
the phonon fluctuations, which holds for small enough ip and is thus justified for an analysis of a continuous transition. 
We also assume replica symmetry in the value of "0 which extremizes the partition function, so that we present below 
an effective free energy of ip which does not involve replicas and which we believe describes the NA transition in a 
random network. 

The integration of v is sketched in the Appendix, and yields two distinct theories, depending on whether or not 
the system exhibits 'soft' elasticity (A = 0). The effective free energies are: 

F^^i ^rol^P - AolVI^ + \Pom' (A = 0) 
V ~ \ '^ta\M' + iPaIM* (A^O) 
which can also be presented as a single crossover expression covering both cases: 

^ = ^r{AM' + A(A, V') + fmm' (4.7) 



The free energies above comprise one of the primary results of this work. In Eq. (4.7) the coefficients t(A) and /3(A) 
are continuous analytic functions of A, while A(A, tp) is a non-analytic function which gives the cubic term for A = 0. 
This term is responsible for the qualitative results we find. In this work we give only the detailed expression for the 
cubic crossover function A(A, ■0), since this term alone determines the qualitative aspects of the transition. 

Soft nematic elastomers (A = 0) — The renormalized coefficient tq has the form fo(T — Tjv^), with the new transition 
temperature shifted by the effects of rubber elasticity and random pinning ( |3.9[ ). Besides this obvious effect, the new 
feature of this soft regime is the cubic term, which restores the HLM effect. This term is given by 



1 (skA 



3/2 



and is precisely the cubic term given by Halperin-Lubensky-Mal. Hence, while the coupling of elasticity to the director 
(or gauge) field destroys the HLM effect, the effect is restored if the coupling preserves the set of Goldstone modes 
present in the nematic state, an intuitively pleasing situation. 

Conventional nematic elastomers (A 7^ 0) — In a network formed in the nematic state there are no soft elastic 
phonons and no Goldstone modes, and the effective Landau-de Gennes free energy contains renormalized quadratic 
and quartic coefficients. There are numerous contributions to the quartic terms, but the qualitative nature of the 
renormalized quartic follows from examining the behavior of the non-analytic crossover function, given by 



A(A,V')= ^ 



127rK3/2 
1 

127r/^3/2 



AAbi + (61 + 62 



2\ , 3 2 4^1 +^2 \W 



61 Ai/2 



(4.9) 



where the first inequality defines the validity of expanding the 3/2 power, and the second inequality allows us to 
simplify the results for a nearly-soft system (A ~ 0). 

Hence, the non- analytic term yields a negative quartic coefficient when expanded for small |0| at non-zero A. For 
small enough A this term becomes very large and the expansion is only relevant for extremely small \ip\. At the 
point where the expansion becomes invalid the quartic term must be replaced by a cubic term, and a tricritical point 
results. Alternatively, one may examine the free energy above, Eq. ( |4.7D , to find a tricritical point at A, given (for 
small A) by 

where P is the bare quartic coefficient. The NA transition is of second (first) order for positive (negative) /3a- 

Thus, for small enough A the transition is expected to be of first order. This correction is most important in those 
systems which are already close to a tricritical point, /3 > [typically those with a relatively small nematic range, 
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Tjjv — Tj^A, where Tjjv is the isotropic-nematic transition temperaturecafl] . We have checked numerically that, for a 
wide range of realistic elastic constants, the correction to f3 ranges up tp. 0.1 kgTq'^ ~ /3, where qc is the microscopic 
(large wavenumber) cutoff that defines the coarse-graining of the theoryQ. As noted in the Appendix, the only effect 
of disorder is to slightly increase the transition temperature, which is sensible because crosslink sites locally enhance 
the smectic order. Our effective free energy assumes a replica-symmetric ground state. In the Appendix we note that 
the replica-symmetric solution is stable for \tlj\ — 0, while in the smectic state there is a crossover at small \tp\ to a 
replica-non-symmetric state, which in this context refers to the effect of the smectic order on the background phonon 
spectrum. An analysis of the further effects of disorder on the low temperature state is beyond this scope of this 
work. 



V. CONCLUSION 



A. Experimental consequences 



We conclude by first discussing some experimental consequences and signatures of this work. The primary qual- 
itative prediction concerns the order of the transition. The preparation of networks of varying degrees of 'softness' 
(which may be controlled by, for example, changing the degree of order which is frozen into the smectic state) allow 
for tests of the predicted tricritical behavior between first- and second-order nematic-smectic-A transitions, which 
may be probed by, for example, specific heat experiments. 

It would also be interesting to test the crossover to a 'glassy' phase characterized by replica-symmetry breaking. 
X-ray scattering to probe the smectic density wave should yield this information, in the form of a crossover from 
a Landau-Peierls' form to a structureJactor with temperature-indepeadent logarithmic correlations, similar to the 
vortex glass considered by Korshuno\i23 and Giamarchi and Le DoussaS. 

A cornerstone of our treatment involves the coupling between the smectic phase variable and the elastic strain, in 
Eq. (1.2), which we have asserted to be a dynamic effect. A signature of this would be a characteristic timescale 
T* associated with hopping of mesogens oyer the smectic barrier. In strong smectic phases this timescale should be 
essentially infinite so that we may use Eq. (1.2), but it should be present in weak smectics and emerge in measurements 
of the complex rheological response function. 



B. Critique and outlook 

In this work we have examined the onset of the nematic-smectic-A transition in a smectic elastomer, and delineated 



two types of behavior (summarized in the free energies, Eqs. 4.6). Soft nematic elastomers should display a first order 
transition, due to the HLM effect induced by the spectrum of Goldstone modes which are a combination of director 
and strain fluctuations. Conventional 'hard' elastomers, on the other hand, should have a continuous transition, with 
a tricritical point as the soft limit is approached. 

In performing our calculations we have assumed the system is in the type-I limit, so that a mean field treatment 
is valid. Since director fiuctuations have an additional mass due to the coupling to the elastic network, we believe 
this limit is safe for smectic elastomers. We have treated disorder within the replica formalism and found that, at the 
level of the onset of the transition, the only effect is to stabilize the smectic state. A preliminary analysis suggests 
that the effects of disorder are certainly important deep in the smectic state. 

Hence a natural extension of this work concerns the nature of the low temperature phase. The effects of disorder 
are two-fold: (1) Translational symmetry is broken and the smectic phase is locally pinned to the disorder, which 
should destroy the Landau-Peierls transition in favor of true long-range order of the smectic density wave; and (2) 
disorder can be strong enough to destroy the long-range order itself. Which, if either, of these effects wins out is 
an interesting question, analogous to the effects of disorder on other systems with continuous symmetry, such as the 
XY-model or flux lattices in superconductors. 

A second, more speculative, result of our work is the suggestion that the coupl ing proposed previously to describe 



the energy penalty for sliding the smectic phase relative to the crosslinks, Eq. (1.2), is actually present only when 
one considers sufficiently short timescales (which, practically, may still be of order of years). This leads to predictions 
of rheological effects, as well as possible glassy phases to 'freeze' this term in, which we leave for the future. This 
interesting suggestion is based on the notion that the smectic layers are actually 'phantom' layers, and their motion 
need not correspond with center-of-mass motion of the mesogens. 

We are grateful to D. Khmelnitskii, T. Lubensky, and M. Warner for helpful conversations. This work has been 
supported by Unilever-PLC (EMT) and the EPSRC (PDO). 
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APPENDIX A: INTEGRATION OF FLUCTUATION MODES 



In this Appendix we sketch the steps to integrat e th e network phonon field v(r) from Eq. ([1.3|). RecaU first our 
coordinate system v = {vz,v±,vt), as given in Eq. (O). We begin with the energy governing the elastic degrees of 
freedom, Eq. ([f.5|), after integrating out the director degrees of freedom: 



^ = ^E / v'^(q).M(g).v''(-q)-r5] [d\\r\\^''\cosqoz-{v"-v'), 



(Al) 



where the matrix M (q, is given by 

/ Czzqlgniq) + Cgi^i + Kaiq^q\ qzq±{C^z + Ka^q^) 





1 



V 





C^^qlgniq) + C'^^ql + Ka^q^ql 

Ciq\gn{q) + C'^2ql + Ka2q^ql 



(A2) 



where gn{q) = + Dq(iP) + Kq^ and the new renormalized elastic moduli are 

Czz 



C3 + fC'l + |C2 + f C4 



Ca^z = (5i+i5o)a± + (6?-fol)/4 ai = A + (61 - 62) V(4foi) 
a2 = A + (61 + 62)V(46i) = C3 - fCi + iC72 - §^4 + C5 

C^i = Afci + Doai C^2 = A6i + i:ioa2. 

Recall that A = C5 — 62/(4^1) is the deviation from an ideal soft nematic elastomeril@. 

We first rescale v, pulling out the determinant of M(q), which leaves the following replica partition function: 



(A3) 



(A4) 



where 



a ^ I a 1 a,b 



qo 



(A5) 



Note that the matrix ]V[(q) becomes a differential operator in real space, with the square root defined with reference 
to the operator in Fourier space. r-irn 

To calculate F^^ (ip) we use the Gaussian Variational Method (GVM)EZlO. Since this method is Ji2ilt|kiiown and is 
essentially the Hartree approximation, we sketch the procedure and refer the reader to reference£3'LHlEil for further 
discussion. We begin by assuming a trial Hamiltonian, 



^* = ?E / w'^(q).G-j,(q).w''(-q). 



a,b ■'q 

The free energy for the w integration satisfies the inequality 

Fw < Fvar = Ft + {H ~ 'Htjt, 



(A6) 



(A7) 



where the average is taken with the trial Hamiltonian Tit- The variational free energy Fyar, which must then be 
minimized over the matrix Ga.b(q), is 



Fva 



^(TrlogG-TrG) + n -Pg J d\\r\\i^''\exp (^-^q^k^T J^BaA^l)^ 



(A8) 
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where 



Ba,h(q) = Zm(q)-[G„.„(q) + Gb,f,(q)-2Ga,h(q)]-z„,(q). (A9) 

Here we have defined z„i{q) = M^^/^(q) •z(q). 

Next we calculcite the self-consistency conditions 5Fu(^ij- j SGa.h — 0- For ci replicci-syninietric Rnscitz, G^^^ — clq^ a,b ~^ 
ai, the self-consistency conditions yield = 1 ^nd 

ai = -2z„,(q)z,„(q)r7<7o'yrf'HV''l (AlO) 
7 = exp \-ql k,T f M-,\q)\ . (All) 



In obtaining this relation we have assumed that the the smectic instabihty is satisfactorily described by a replica- 
symmetric ansatz jV'^l = \ip\- Henceforth we discard free energies terms of higher than linear order than n, which 
vanish in the n — )■ limit. Moreover, we expect that in the low temperature phase (|V'''I)p is replica-symmetric, and 
glassy behavior manifests itself in broken replic£usjimmetry of the phase variable [i.e. the layer spacing u(r)] as in, 
for example, the XF-model in a random ficldEjE3'E2l. 

Next we identify the free energy using Eq. (B.7), and examine the stability of the nematic phase against uniform 



smectic order. Hence we examine the stability of = fcBTlim„^o(l ^ (-2") p)/n. The final result is: 

^ - irlV'P + + ksT log5„(9) + ifc.T^logdet M (q, ^)-T jV^'l (1 - 7). (A12) 

There are four g-integrals to perform: one (/i ) from log g„ (q) , which arises from director fluctuations and from the 
logarithm of the prefactor of the determinant of M. These are identical, and exhibit a partial cancellation [a factor 
of 1 from the director fluctuations and a factor of —3/2 from the determinant of M]; two more integrals {I2 and I^) 
emerge from logdetM, and one from {I4) Mj^,^, which enters in the exponential in 7. 

Ii may be performed exactly and expanded to yield even powers of ip. We take I2 to be the integral of the logarithm 
of the single diagonal of M involving C4 (that is, Mft). This has a contribution non-analytic in ^jj for A ^ 0, since 
in this case C52 ~ l-iAP and the integrand in the q-integral is singular as — > 0. This yields the cubic term A(A, if)), 
as well as many terms even in tp. is then the logarithm of the determinant of the remaining 2x2 sector of M. 
This is also non-analytic in ^/^ at A = 0, but yields terms of order jV'l^, "0* log IV'li which we ignore compared to 
the effect of tj}^ contribution (these terms can be shown to lead to subdominant terms for small A). 

The last integral, appearing in 7 = e^*, must be handled with care. We are interested in expanding this integral 
for small ip for both the soft (A = 0) and hard (A ^ 0) cases. For A,'(/; — > the integrand is singular as 0. 
This singularity yields a logarithmic contribution to the integral, so we write 



where /J is a constant and fo{ip) is an analytic function of ip. Hence the term involving 7 in Eq. ( A12| ) contributes a 



term in the energy psoportional to -0^^^, multiplied by an exponentially small prefactor. We find 77 to be essentially 
the Caille exponentllll, r] = q^k gT ^Jbi/K / [%ttC) ^ where C is a combination of rubber elastic constants. Since 77 may 
be of order 0.5 — 2.0, we ignore this term compared to the cubic term. If 77 is sufficiently small we may include this 
term in a Landau expansion, but we expect the exponentially small prefactor to render it negligible. In the case where 
A is non-zero the integral is analytic and yields even powers of '0. Therefore, the only qualitative effect of disorder is 
to increase the transition temperature. 

The final step of the calculation is to ensure stability of the Keelica-symmetric ansatz for the v integral. This is 
done by examining the replicon mode, identified for the GVM asEll 

A = 1 - r|V^|2g4fc^T7 / (M-i)2, (A14) 



where A is the eigenvalue of the most unstable fluctuation mode about the replica-symmetric solution given by (ao, ai). 
We note that this is actually the eigenvalue of the fluctuation mode of the 7V„ ^ 00 limit version of the theory, where 
Ny is the number of 'color' cornoonents of the field v, in which limit the GVM (or Hartree) approximation becomes 
the exact saddle point integraEil. While A is calculated in this limit, it may or may not apply to the physical case 
= 1, but hopefully yields intuition about the correct behavior. 
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Since A = 1 in the high-temperature phase (because = 0), the rephca-symmetric solution is stable as far as 
understanding the onset of the smectic transition. Analysis reveals an L-divergence in J^{'M.~^)'^. which indicates 

a critical temperature (in the smectic state) at which the disorder is relevant and we must consider replica-non- 

symmctric states. Here L is the system dimension. Ignoring numerical factors, this 'glass' transition temperature T* 
is given by the temperature very slightly below the NA transition at which the order parameter attains the value -0*, 
given by 

\r\- -oi^] , (A15) 



where C ~ kgTNx is a combination of rubber elastic constants. Since tp* is inversely proportional to the system 
dimension L, the window within which the replica-symmetric solution holds is expected to be very small. Thus, for 
understanding the properties of realistic systems, one must examine the properties of the disordered low temperature 
state within the framework of replica-symmetry-breaking. 
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